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ABSTRACT 

We present late-time Hubble Space Telescope (HST) imaging of the fields of six Swift gamma-ray 
bursts (GRBs) lying at 5.0 < z < 9.5. Our data include very deep observations of the field of the 
most distant spectroscopically confirmed burst, GRB 090423, at z = 8.2. Using the precise positions 
afforded by their afterglows we can place stringent limits on the luminosities of their host galaxies. 
In one case, that of GRB 060522 at z = 5.11, there is a marginal excess of flux close to the GRB 
position which may be a detection of a host at a magnitude Jab ~ 28.5. None of the others are 
significantly detected meaning that all the hosts lie below L* at their respective redshifts, with star 
formation rates SFR < 4 M0yr _1 in all cases. Indeed, stacking the five fields with WFC3-IR data 
we conclude a mean SFR < 0.17 M Q yr _1 per galaxy. These results support the proposition that the 
bulk of star formation, and hence integrated UV luminosity, at high redshifts arises in galaxies below 
the detection limits of deep-field observations. Making the reasonable assumption that GRB rate is 
proportional to UV luminosity at early times allows us to compare our limits with expectations based 
on galaxy luminosity functions derived from the Hubble Ultra-Deep Field (HUDF) and other deep 
fields. We infer that a luminosity function which is evolving rapidly towards steeper faint-end slope 
(a) and decreasing characteristic luminosity (L*), as suggested by some other studies, is consistent 
with our observations, whereas a non-evolving LF shape is ruled out at > 90% confidence. Although 
it is not yet possible to make stronger statements, in the future, with larger samples and a fuller 
understanding of the conditions required for GRB production, studies like this hold great potential 
for probing the nature of star formation, the shape of the galaxy luminosity function, and the supply 
of ionizing photons in the early universe. 



1. INTRODUCTION 

The drive to locate and characterize the earliest grav- 
itationally bound objects in the Universe - the first 
proto-galaxies, population III & II stars and the black 
holes they produced - is a central theme of contem- 
porary cosmology. It is intimately tied to the quest 
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to understand how the Universe at large was reionized, 
quite possibly by the ultra-violet (UV) flux produced by 
these first, and still enigmatic objects. Major invest- 
ments of time on large facilities continue to be spent 
on these ambitious goals, but progress has been hard- 
won. Historically, quasars were the key high-redshift bea- 
cons, thanks to their great luminosities, but the space 
density of bright quasars drops above z ~ 4 and to 
date the most distan t quasar is located at z = 7.1 
(|Mortlock et al.l |2011[ ). Galaxies are, of course, much 
more common than quasars, and high-z galaxies are 
likely to be present in moderate numbers in very deep 
near-infrared (nIR) observations. Indeed, recent nIR 
observations in the Hubble Ultradeep Field (HUDF) 
have unveiled > 100 candidate z > 7 galaxies (e.g., 
Bouwens et"aHl2011d : [Bunker et alj|2010t iMcLure et al.1 
20101 ) discovered via the Lyman-dropout technique. The 
challenge of studying these sources lies in their ex- 
treme faintness (often >28th magnitude), which gen- 
erally means only photometrically-derived redshift esti- 
mates are possible. Furthermore, distinguishing young, 
high-z galaxies from old or dusty populations at moder- 
ate redshifts, or even Galactic brown dwarfs, becomes 
increasingly difficult when approaching the sensitivity 
limit of the data. Although spectroscopic redshifts have 
been determined for a small number of z > 7 galaxies 
(|Vanzella et al.ll2011h lOno et al.ll2012[ ). this is only pos- 
sible if the galaxy is a strong Lya emitter (and the emis- 
sion is not completely absorbed by a partially neutral 
intergalactic medium), and requires a major investment 
of time for each target. Even in the era of extremely 



2 



large telescopes, many such sources may be too faint for 
direct spectroscopic redshift confirmation. 

The extreme luminosity of gamma-ray bursts (GRBs) 
makes them potentially powerful probes of the early Uni- 
verse, a utility which has been wide ly touted since the 
discovery of the first afterglows (e.g.. iWiiers et al.l[T998t 
lLamb &: Reichartll2000HTanvir fc Jakobssonll2007t ). Cur- 
rent technology has the capability to detect the prompt, 
and afterglow emission for bright G RBs out to z ~ 
20, should GR Bs e xist at this ep och (|Gou et al.1 120041 : 
iRacusin etHI 120081: iBloom et alfl2009f ). and they have 
several advantages over alternative methods for the de- 
tection and study of high— z objects. In the first in- 
stance, the intrinsically smooth power-law spectra of 
the afterglows makes them ideal backlights for absorp- 
tion diagnostics: not only providing redshifts, but in 
principle chemical enrichment, hydrogen column densi- 
ties (e.g.. iVreeswiik et~aT1 120041: lProchaska et all 120071 : 
iFvnbo et al.ll2009l). extinction and dust laws in the hosts 
(e.g.. IZafar et a" I [20llat ISchadv et all l2012h . and even 
probi ng the state of the i ntergalactic medium (IGM) 
(e.g. lMiralda-EscuddH99l iMcQuinn et all 120081 ). Sec- 
ondly, they pinpoint the positions of their hosts, and 
once the afterglow has faded we can search for the host 
galaxy with a redshift in hand from spectroscopic (or 
photometric) observations of the afterglow. Finally, long- 
durat ion GRBs are produced during core-col l apse events 
(e g IHiorth et al.l[2003al iStanek et alll2003t iPian et al.1 
12006( 1. and so trace the locations of at least some com- 
ponent of massive star formation. Importantly, GRBs 
allow absorption redshift measurements for galaxies at 
times too faint to be seen in deep im aging even with the 
Hubble Space Teles c ope (HST) (e.g., iBerger et all [20021: 
IHiorth et al.H2003bt iThone et al.ll201fL 

The early populations of massive stars and the proto- 
galaxies in which they reside are thought to be ma- 
jor producers of ionizing photons, and also give rise to 
GRBs. Hence, mapping GRB space density as a func- 
tion of cosmic time should trace this early star forma- 
tion and so provide a key to assessing its contribution 
to the reionization of the IGM at z > 6. Since GRB 
progenitors are individual stellar systems, their hosts 
should sample the whole (star-forming) galaxy luminos- 
ity function (LF), rather than just the bright end, avoid- 
ing the limitation inh erent in flux limited samples (e.g., 
Uakobsson et al.ll2005l ). This is particularly important at 
very high- redshift, when galaxies were small and faint, 
and only the tip of the LF can be probed directly. 

However, GRB-based studies are subject to their own 
difficulties. As with any sources beyond z ~ 6, in ad- 
dition to the increasing luminosity distance, we have 
to contend with the difficulties of working in the nIR. 
In addition, because GRB afterglows fade there is only 
a narrow window of opportunity in which observations 
can be pursued, so signal-to-noise cannot be built up 
over many nights of observation. Hence, a major con- 
straint in their detection and study at high-z is the 
availability of large aperture telescopes which can re- 
spond rapidly with appropriate nIR instrumentation. It 
is also clear that high-redshift GRBs detectable to the 
limits of the Swift satellite are rather rare, and there- 
fore many burst s have to be observed to find the few 
at z > 5 (e.g.. Uakobsson et al.l l2012h . Finally, stud- 
ies of GRBs at lower redshift have indicated that they 



are generally found in small-to-moderate size, low-to- 
moderate metallicity galaxies with high specific star- 
formation rates (= SFR/M^w e. g.. IChristensen et al.1 
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2004t iCastro Ceron et all 120061: iModiaz et al.l . _ 
Chen et al.l 120091 : iSavaglio et al.1 120091 : iSvensson et al.1 
2010ft . and ra rely in intensively star-forming far-IR- 



bright galaxies dTanvir et al.l l2004t [ Le Floc'h et aLl l2006l: 



iMichalowski et al" _ 20081). Coupled with theoretical ar- 
gume nts (e.g., lYoon fc Langerl l2005t iWooslev &: Hegerl 
|2006| ) and spectroscopi c studies of their env ironments 
(e.g.. lFvnbo et afll2009l : ILevesque et al.ll2010al ). this has 
led to the suggestion that GRBs may be more com- 
mon and/or brighter at rela tively low metallicities (e.g., 
Wolf & Podsiadlowski 2007). Such a dependence could 
potentially help explain the apparent increase in GRB- 
rate to SFR ratio between z ~ and z ~ 5 seen in 
various studies (e.g..lNataraian et al.ll2005l : lYiiksel et al.1 
120081: iRobertson fc Ellis! 12012ft . However, the picture is 
not a simple one: there is a selection bias against finding 
heavily extinguishe d afterglows (e.g.. iPerlev et al.l l2009t 
iKruhler et alll20i"ll) . and at least some GRBs appear to 
be produc ed in very dusty, massive star-bursting galax- 



ies (e.g.. iChen et al.l 
Hashimoto_et_a!Jl20I0; 



2010; IKupcu Yoldas et all I20T0I: 



Svensson et al.ll2012f ) . whilst oth- 



ers have been foun d in solar or even super-solar metallic- 
ity en vironments ()Levesaue et al.l l2010bl : ISavaglio et al.l 
1201 If) . Furthermore, comparison of GRB luminosity and 
host metallicity at z < 1 reveals no obvious correlation 
(jLevesque et al.lfeOlOd ). In any event, it is likely that 
galaxies with properties similar to typical GRB hosts pre- 
dominate at early times, and indeed the small number of 
higher-redshift afterglows for which accurate metallici- 
ties have been derived show a range that matches well 
predictions from galaxy evolution simulations that in - 
clude no bias in GRB production (jPontzen et al.ll2010f ). 
Thus, we expect GRBs to be a good tracer of the bulk 
of high redshift star f ormation (|Fvnbo et al.ll2006l 120081: 
iKocevski et al.l [20091 ) . and we proceed under that as- 
sumption in this paper. 

Several very high-z bursts have been identified by 
Swift to date, and have begun to realise their poten- 
tial as probes of the early universe. The afterglow of 
GRB 050904 at z = 6.3 was brighter th an magnitude 
J = 17.5 even a few hours after the burst (jHaislip et al.l 
2006), and spectroscopy would have been routine with 
a 4-8 m class telescope. Indeed, ultimately the red- 
shift measurement came from a Subaru spectrum taken 
some 3 days after the event ([Kawai et al.l l2006h . Even 
the intrin sically rather faint afterglows of GRB 080913 
(i = 20: IGreiner et al.l [2001) or GRB 090423 (J = 21; 
iTanvir et al.l 20091) allowed for spectroscopic redshifts 
with 8 m instrumentation, despite the data being ac- 
quired more than 12 hours after the burst in the latter 
case. 

Until now, the most distant GRB host detected 
(with Spitzer) was that of GRB060510B at 2 = 4.94 
(jCharv et al.ll2007T ). In this paper we present deep, late- 
time Hubble Space Telescope (HST) observations of the 
fields of six of the most distant gamma-ray bursts, all at 
z > 5. Only one host galaxy, that of GRB 060522, could 
be marginally detected, and the deep limits for the others 
imply that all lie close to or below the characteristic lu- 
minosity, L* at their respective redshifts. As we show in 
this paper, the non-detections of the hosts place impor- 
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tant constraints on their total luminosities and star for- 
mation rates. Since long-duration GRBs trace (at least) 
some component of star formation, a survey of even a 
small number of GRBs at high redshifts provides a cen- 
sus of the locations and galactic environments of star 
formation at early times. This means that GRBs have 
the potential ultimately to constrain the faint end of the 
galaxy LF at z > 7, which is crucial for understanding 
the reionization of the Univ erse, thought to occu r pre- 
dominantly at 8 < z < 12 (jKomatsu et al.l 120111 ). UV 
photons produced by massive stars are widely c onsid- 
ered t he most likely driving force for reionization (jLoebl 
2009; Fa ucher-Giguere et al.l [2009h . but even ultra-deep 
surveys cannot currently quantify the (likely dominant) 
contribution of intrinsically faint galaxies. 

Throughout this paper we use the AB-magnitude 
system, and adopt a A-CDM cosmology with Ho = 
72kms- 1 Mpc- 1 and Q M = 0.27, Q A = 0.73. 

2. GRBS IN THE SAMPLE 

Our sample comprises six of the most distant bursts 
detected by Swift, including all five bursts detected to 
date that have firm spectroscopic redshift s above z = 

5. They are GRBs 050904 (z = 6.29; iKawai et al.l 

l200ll . 060522 (z = 5.11: ICenko et all 120061) . 060927 
{z = 5.47: iRuiz-Velasco et al.1 12007D 0809 13 (z = 
6.73; IGreiner et al.l 120091: IPatel et al.l [2010h. 090423 



2009) 

20TTF . 



(z = 8.23; iTanvir et al.l 120091 iSalvaterra et al. 
and GRB090429B (z phot = 9.4; iCucchiara e~ 
Here, we briefly summarise details of each GRB, and 
the available HST observations of each field, which were 
mostly obtained after the afterglow should have faded be- 
yond detectability (the possible exception is GRB 050904 
as discussed below). The majority of our data comes 
from the new Wide Field Camera 3 Infrared channel 
(WFC3 IR). In all of these cases we used the standard 
flat-fielded observations from the ilST archiveQ, and cor- 
rected these for geometric distortion, creating a stack of 
individual images using multidrizzle. The final pixel 
scale for these well dithered observations is set as 0'.'05 
per pixel, with the pixf rac set to 0.7. For observations 
with ACS we drizzled to the same final pixel scale, which 
broadly retains the native pixel scale of the instrument. 
For GRB 060927, we also present NICMOS observations 
of the host, as detailed below. Cut-outs of the images 
around the locations of the GRB positions are shown in 
Fig.ffl 

The exact locations of each burst on the HST images 
were found relative to field sources using ground-based 
images of their afterglows. Details of the data obtained 
for each burst and the rms accuracy of the astrometric 
calibration are summarised i n Table [U 

For photometry, following iBouwens et al . (2010aj), we 
employed O'M diameter apertures centred at the loca- 
tions of the afterglows (for reference, 0'.'4 corresponds to 
a physical scale of « 2.1 kpc at z = 7). We note this is 
appropriate since GRBs are generally fo und to be close 
to th e UV-brightest regions of their hosts (|Fruchter et al.l 
Hxl), and the sizes of high-z galaxies are found to be in- 
trinsically small (e.g. 0.7 ±0.3kpc was found for z « 7-8 
galaxies in the HUDF by lOesch et alll2010f ). The prior 
knowledge of the exact locations of the GRBs allows us to 

14 archive.stsci.edu 



adopt a 2cr excess flux as a reasonable detection thresh- 
old, which is unlike blind searches where a much higher 
level of significance is required to both ensure confidence 
that the source is real and also to reliably constrain the 
photometric redshift. For comparison to galaxy samples, 
and GRB hosts at other redshifts, it is relevant to con- 
sider the appropriate aperture correction^, and these 
are included in the limits shown in Table [2] Our errors 
are determined in each case from the variance within a set 
of apertures of equal size to our source aperture, placed 
at random blank sky locations in the field surrounding 
the source. 

The results of the photometric analysis, and inferred 
limits on host properties, are given in Table [2j Note 
that we report the measured sky-subtracted flux densities 
even when they are not significant detections, or indeed 
are formally negative. This is to allow comparison with 
models, as is done in Section 3. 

2.1. GRB 060522 

Spectroscopy of GRB 060522 was performed by Keck 
roughly 12.5 hours after the burst, and revealed a strong 
break ~ 7425 A in the spectrum, interpreted as Lya at 
z = 5.11±0.01 ()Cenko et al.ll2006[ K A search for the host 
galaxy with the Spitzer Space Telescope found no detec- 
tion to flux densities of 0.2 fj,Jy at 3.6 /im and 2.4 /xJy at 
5.8 ^m, implying that its r est-frame optical luminosity 
lies substantially below L* (|Charv et al.ll2007f ). 

We obtained observations of GRB 060522 with 
WFC3/IR using the F110W filter. To ascertain the pre- 
cise location of the burst on our HST image we performed 
relative astrometry between our WFC3 observations and 
an image obtained at th e Telescopio Nazional e Galileo 
(TNG) on 22 May 2006 (jD'Avanzo et al.ll200ll . At the 
location of the afterglow we measure a flux density of 
7 ± 4nJy, corresponding to a 2a (aperture corrected) 
magnitude limit of FllOW(AB) > 28.13. Although for- 
mally this is a non-detection at 2a it is possible that 
some host flux is contributing within the aperture. In- 
deed visually it appears there is a somewhat more signifi- 
cant excess of flux which is offset slightly south-west from 
the GRB location (by about 0'.'3, or ~ 2 kpc at z = 5.11), 
and if we increase the aperture size to 0'.'6 diameter then 
the flux density becomes 12±5nJy. This would represent 
a marginal detection, and corresponds to an apparent 
magnitude of F110W(AB)w28.5. While this is a plausi- 
ble magnitude for a host (corresponding to Migoo(AB) « 
— 18), the probability of a chance alignment with an un- 
related object at such faint magnitude levels is also non- 
negligible. Specifically, cons idering the F110W num ber 
galaxy number counts from iThompson et alj (|1999f ) we 
estimate « 5% of random locations will have a galaxy 
of magnitude F110W(AB)< 28.5 within 0'.'5. Thus in a 
sample of six, as we have, there is a ~ 25% chance of 

15 i.e. the limits in this table are given by Mr Im = 
23.9 — 2.5 log (flux + 2<t) + &p corr . The aperture cor- 

rections are calculated from the encircled energy (e.g. 
http://www.stsci.edu/hst/wfc3/documents/handbooks/currentIHB/ 
cUYjrUY.htmlJ an3 nave values for W-b'U3/lK 51 
apios = 0.31,apno = 0.33, api25 = .39, api60 = 
0.52 for ACS/WFC ap 85 o = 0.27 l lSirianni et al.l 
20051) ■ and for NICMOS /NIC 3 apim - 0.60. 



http:/ /www. stsci.edu/hst/nicmos/documents/handbooks/ cur- 
rent JNhJW/c04_imagmg.t>. 5. html 
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TABLE 1 

Summary of burst sample and astrometric ties 



GRB 


z 


z-ref 


Obs. summary 


Astrometry 


Accuracy (") 


060522 


5.11 


1 


WFC3:F110W 


TNG 


0.06 


060927 


5.47 


2 


NIC3:F160W, WFC3:F110W 


VLT/FORS2 


0.07 


050904 


6.29 


3 


ACS:F850LP 


Gemini-S/GMOS 


0.06 


080913 


6.73 


4 


WFC3:F160W 


VLT/FORS2 


0.08 


090423 


8.23 


5 


WFC3:F125W,F160W 


VLT/HAWK-I 


0.03 


090429B 


9.4* 


6 


WFC3:F105W,F160W 


Gemini/NIRI 


0.06 





Note. Referen ces { 1) ICenko et al.l BOOBI). (2) IRui z-Vclasco ct al. 1200711. (3] 

IKawai et all $200® ). (4) IPatel et al.l H20101 ). (5) ITanvir et al.l H20091 ). f6) ICucchiara et al 
(2011). *Photometric redshift estimate. 



F850LP 060522 




F110W 060927 



060927 F110W 




080913 



F160W 090423 F125W+F160W 090429B 



F160 





Fig. 1. — A mosaic of the host galaxy fields of the z > 5 GRBs in our study. The filters of the observations are indicated in each panel, 
while the cross-hairs represent the location of each burst as determined from astrometry with ground based images containing the afterglow 
(see Table [T] for more details). The images have been lightly smoothed to allow the eye to see fainter features. 



an unrelated galaxy being close enough to be possibly 
mistaken for a host in at least one case. For the sake of 
consistency with the other bursts we will continue to use 
the Of.' 4 aperture limit in the remainder of the paper. 

2.2. GRB 060927 

GRB 060927 was detected by Swift, and its optical af- 
terglow was initially found in rapid, but unfiltered Rapid 
Optical Transient Source Experiment (ROTSE) obser- 
vations. An extensive f ollow-up campaign (reported in 
iRuiz-Velasco et all I2007D revealed the source to be an 
i?-band drop-out, well detected in the i-band. A sub- 



sequent spectrum obtained at the Very Large Telescope 
(VLT) showed a faint continuum redward of ~ 8000 A, 
as well as a weak Si II absorption feature, consistent with 
a redshift of z = 5.47. 

Our ffST/NICMOS observations of the field of 
GRB 060927 were taken using the F 160W filter 
and were reduced as described in iFvnbo et al.1 
(2005). GRB 060927 was additionally re-observed 
with WFC3/IR using the F110W filter. Astrometry 
was performed via observations made with the VLT on 
30 September 2006. To establish the position on the 
NICMOS images (which have a narrower field of view 
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than WFC3) we opted to perform relative astrometry 
directly to the WFC3 images. The resulting rms 
astrometric accuracy is <(y.'01, and so does not impact 
the overall error in the astrometric solution described 
above. 

At the location of the afterglow the photometry 
in each of the NIGMOS and WFC3/IR observations 
yields limiting magnitudes of F160W(AB) > 27.75 and 
FllOW(AB) > 28.57 respectively. 

2.3. GRB 050904 

GRB 050904 was the first GRB at z > 6 to be lo- 
cated. Its optical/ IR afterglow was initially found by 
lHaislip et al.l (|2006l ) . who derived a photometric redshift 
of z = 6.39 based on the strong spectral break between 
the and i- and z-bands, coupled with nIR observations 
sho wing a blue spectral slo pe redward of the break (see 
also iTagliaferri et al.ll2005l ). The afterglow was intrinsi- 
cally extremely bright, amongst the brightest observed 
for any burst, making late time spectroscopy feasible 
with Subaru. This provided a measurement of the ab- 
sorption redshift, hydrogen column density, and meta l- 
licity of the host (jKawai et alJl2006t [Totani et al.ll2006f ). 
It remains the most distant burst for which all of these 
diagnostics are available. 

A search for the host galaxy of GRB 050904 was con- 



ducted by iBerger et aLT (|2007D using both HST and 
Spitzer, which placed deep limits on any host emission. 
The first HST epoch, at which time both ACS/F850LP 
and NICMOS/F160W images were obtained, was carried 
out only about three weeks post-burst. No significant 
flux was detected in F850LP, and the faint detection in 
the F160W image was shown to be due to residual after- 
glow c ontamination since it was absent in a later F160W 
epoch (jBerger et al.l 120071 ) . 

For our analysis, we re-reduced the ACS data for 
GRB 050904 (since we are primarily interested in the 
rest-frame UV luminosity close to Lya, this filter is the 
more relevant). Astrometric tying of the afterglow to 
field sources was done utilizing a z-band image obtained 
from Gemini-South on 7 September 2005. At the loca- 
tion of the afterglow in the F850LP frame we measure a 
flux density of — 5 ± I7nJy, corresponding to a 2a limit 
of F850LP(AB) > 27.50 (or 27.06 at 3a) This is in 
good a greement with the limits reported by IBerger et al.l 
(|2007| ). However, since the Lya break lies within the fil- 
ter bandpass at this redshift, we must account for flux 
lost due to IGM absorption (i.e., the effective filter width 
is narrower for this host), and so we conclude a corrected 
magnitude limit of F850LP(AB) > 26.86. 

2.4. GRB 080913 

GRB 080913 was identified as a high redshift candi- 
date based on photometric observations wi th GROND 
show ing the burst to be an z-band drop out (|Rossi et al.l 
2008). Deep, red spectroscopy from the VLT showed 
a strong spectral br eak, interpreted as Lya at z ~ 6.7 
(jGreiner et al.ll2009fl. A furthe r detailed analysis of the 
spectrum bv lPatel et al.l (|2010f ) revealed a single absorp- 
tion line of Sill at z = 6.733, and we adopt this as the 
redshift of GRB 080913. 

We obtained HST observations with WFC3/IR in the 
F160W filter. To tie the astrometry of our HST obser- 
vations we utilized images obtained from FORS2 at the 



VLT on 13 September 2008 in the z-band, and images 
taken with NIRI on Gemini-North on 14 September 2008. 
We used two images independently to confirm the pre- 
cise location of GRB 080913 on our HST images. For the 
FORS2 images we identified 9 compact sources in com- 
mon, while only 6 sources were usable from the Gemini 
observations. Although the afterglow detections are of 
low S/N, the errors on their centroids are small in com- 
parison to the errors derived from the fit and hence we 
are able place the afterglow to a rms accuracy of 0'.'08. 
Again, no host galaxy is visible in our observations to a 
limit of F160W(AB) > 27.92. 

2.5. GRB 090423 

GRB 090423 was first identified as a candidate high- 
rcdshift object based on its aft erglow being a V-ban d 
dropout, which implies z > 7.5 (|Cucchiara et al.| [20091 . 
VLT spectroscopy with both ISAAC and SINFONI al- 
lowed the identification of the Lyman break, despite 
low S/N, est ablishing the redshift z = 8.23 ± 0.07 
(|Tanvir et al.ll2009l ). This value is in excellent agreement 
with the z = 8.1+5^3 determined from a spectrum ob- 
tained at the TNG (|Salvaterra et al.ll2009H . 

We obtained our first HST observations of GRB 090423 
on 24 January 2010. At this stage we acquired 20 orbits 
in each of the F125W and F160W filters on WFC3. Un- 
fortunately 15 of these 20 orbits were substantially im- 
pacted by persistence from earlier observations of bright 
field sources. This affected the sensitivity over wide re- 
gions of the detector, although does allow us to mea- 
sure the flux at the location of the GRB host (which 
was only mildly affected by persistence). Additional ob- 
servations were obtained in October 2010. All of the 
available observations in each filter were co-aligned and 
stacked via multidrizzle. These observations are ex- 
tremely deep, with limiting magnitudes a factor 1.5 x 
deeper than the WFC3 Early Release Observations of the 
Chandra Deep Field South/Great Observatories Origins 
Deep Survey (CDF-S/GOODS) fields, and only another 
factor of 1. 5 x shallower than the WFC3 observations of 
the HUDF (|Bouwens et al.ll2011cD 

The location of the afterglow on our HST images was 
achieved by tying the astrometry to VLT/HAWK-I ob- 
servations that were obtained approximately 17 hours 
post burst. While the afterglow had faded since its first 
discovery from Hawaii, these deep observations yielded 
similar signal-to-noise, with the added advantages of a 
larger field of view, and a greater number of faint sources 
for comparison. 

Within the HST image there is no obvious source at 
the afterglow location, the measured 2a limiting mag- 
nitudes are F125W(AB) > 30.29, F160W(AB) > 28.36. 
Since both filters sample the rest-frame far-UV, we form 
a weighted average of the two results to provide a com- 
bined flux density measure of —0.15 ± 1.7 nJy, which is 
used to derive the limit on the SFR reported in Table El 
Such a small host is consistent with the non-detection 
of mo lecular gas from the GRB location (St anwav et al.l 

2.6. GRB090429B 

GRB 090429B is the only host within our sample that 
does not have a spectroscopically measured redshift. 
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However, the photometric break between the J- and H- 
bands, coupled with the blue spectral slope between H 
and K provides a best-ht photometric redshift of z = 9.4 
and a robust lower li mit to the redshift of z > 6.5 
(|Cucchiara et al.ll2011l ). As we commented above, the 
smooth power-law spectra of GRB afterglows makes pho- 
tometric redshift estimates generally more reliable than 
they are for galaxies, since the ra nge of intrinsic spec - 
tral variation is much less (e.g., iKruhler et al.l 1201 If ). 
We obtained observations in F606W (ACS), F105W, and 
F160W; the non-detection of a host galaxy in these im- 
ages, including the blue filters, provides additional sup- 
port for the high-z origin for this burst, since the hosts of 
GRBs at z < 3 have so-far always been detected in HST 
optical imaging. Here we consider primarily the F160W 
observation, since this is the only filter redward of Lya at 
the best-fit redshift z = 9.4. However, for completeness 
we also report results for the z ks 6.5 lower limit, which 
allows us to use both the F105W and F160W data. 

We ascertained the location of the burst on the HST 
images via relative astrometry between our first epoch 
-ftT-band observations and those obtained with HST. The 
2(7 limiting magnitude at this location is F160W(AB)> 
27.78. This is shallower than for the majority of the 
bursts in our sample, since only a two orbit exposure was 
obtained. However, despite this the image still probes to 
faint limits comparable to the likely characteristic galaxy 
luminosity, L* , at z — 9.4. 

3. ANALYSIS AND DISCUSSION 
3.1. Limits on host properties 
3.1.1. Star formation rates 

At the redshifts of the bursts in question our nIR 
observations probe rest-frame wavelengths roughly in 
the range 1300-2000 A. Limits on the UV luminosity 
at these wavelengths provide direct constraints on the 
host star-formation rates. A potentially important con- 
sideration in the UV is the effect of extinction by 
dust, which could lead to a significant underestimate 
of the true SFR if uncorrected. However, observa- 
tional constraints from GRB afterglows at high red- 
shifts tiGreiner et al.ll2009tlTanvir et alJ200alZafar et al.l 
1201 lbD suggest that extinction corrections are likely to 
be small. Similarly the blue colours of the z ~ 7 can- 
didate galaxies identified in the deep HST fie lds (e.g., 
IBouwens etld] (l2010bf).lFinkelstein et al.l (|2011ft . but see 
also iMcLure et al.l (|2011[ )) also argue for little dust in 
most early star- forming galaxies. We therefore assume 
dust extinction can be ne glecte d. 

Following lMadau et all (jl998fl we estimate the star for- 
mation rates base d on the UV lumino sity at ^1500 A (see 
also discussion in iBunker et alllMol) : 



SFR = 



Ll500,UV 



8 x 10 27 ergs s^Hz 



— MQyr- 1 (1) 



As discussed in Section 12. 1[ there is a marginal detec- 
tion of what may be the host of GRB 060522 at z = 5.11, 
slightly offset from the burst position, but otherwise none 
of the hosts are significantly detected. The inferred 2a 
limits on their star formation rates are given in Table G2 
and in Figure [2] are plotted as a cumulative histogram of 
upper limits, compared to the SFRs for a sample of z ~ 7 



HUDF galaxies from IBouwens et all (|2011cf ). The limits 
span the range 0.4 — 4 M Q yr _1 , and indicate that the to- 
tal star formation rates in these galaxies are modest. For 
comparison, the median SFR for GRB hosts at z < 1 ar e 
found to be around 1-2 M yr _1 ()Svensson et al.| [2010). 

3.1.2. Stacked limits on host galaxy emission 

The observations of six high-z fields provide the op- 
portunity to stack the resulting data in an attempt to 
provide either a combined detection, or composite limit 
on the observed host luminosity. However, this is com- 
plicated by the use of different instrument and filter 
combinations to obtain the imaging. Therefore, we do 
not attempt to include data taken with different instru- 
ments, and consider a stack of only the WFC3 observa- 
tions. Thus, we omit the observations of GRB 050904, 
and the NICMOS observations of GRB 060927. To per- 
form the stack we first re-scaled the individual images 
(taken in different filters) such that the units of each im- 
age can be considered equivalent. We then computed 
the necessary offsets to overlay the burst positions, and 
re-drizzled onto a single output frame with these offsets 
applied. Simple aperture photometry at this location was 
then performed, with our errors based on the scatter in 
background apertures as before. The mean luminosity 
distance is 68500 Mpc, corresponding to a mean redshift 
of z = 6.82. There is no significant excess flux at the 
afterglow location, and formally the measured stacked 
flux density is 1 ± 3 nJy. Hence, the mean flux density of 
each host galaxy is constrained to be 0.2 ± 0.6 nJy, cor- 
responding to a 2a limit of m^B > 30.7, equivalent to an 
absolute magnitude M\it\ +Z \ < —16.2 and SFR <0.17 
Mq yr _1 , at this mean redshift. 

3.2. Constraints on high redshift galaxy luminosity 
functions 

Although our small sample of high-z GRBs does not 
yet allow us to place strong constraints on the galaxy LF, 
we can test whether the limits on the host magnitudes are 
consistent with them having been drawn from the LFs 
suggested by other studies of galaxy populations with 
redshifts between 5 and 10. 

The recent re-observations of the Hubble Ultra-Deep 
Field with WFC3/IR have revealed a population of z- 
band and V-band dropouts, with colours c onsistent with 
galaxies at z > 7 and z > 8 respectively (|Bunker et al.l 
120101 : IMcLure et alJl20l5 IBouwens et alJl2010a[) . A sin- 
gle candidate z ~ 10 galaxy (J-band drop out) has also 
been identified (Bo uwens et al.ll2011a|) . From these sam- 
ples, assuming they are substantially complete and un- 
contaminated, it is possible to make some statements 
about the form and evolution of the galaxy LF from 
5 < z < 10. In particular, these authors fit their data 
in bins of redshift with th e LFs described by a Schechter 
function (|Schechterlll976l ): 



(x)dx = 4>*x 



c dx 



(2) 



where x = L/L* , with L* being the characteristic lumi- 
nosity of the "knee" of the LF. Here a is the power-law 
slope towards faint luminosities and <p* is a normalisa- 
tion factor. Of course, there are no strong observational 
reasons to expect the LF to have thi s form at high-z, b ut 
it is supported by theoretical work (|Trenti et al.l l2010). 
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Fig. 2. — A cumul ative histogram of the star formation rates (based on the 1500 A magnitude) of a sample of z ~ 7 galaxies (from 
IBouwe ns ct al. 2011c, where we have assumed each galaxy lies at z = 7.1) compared with the inferred 2<r limits on the star formation rates 
from our observations of GRB hosts. As can be seen, the limits attained in several individual cases correspond to galaxies in the faintest 
third of the high-z galaxy distribution (i.e. limits in the case of all our WFC3 observed hosts). This confirms that much high— z star 
formation is taking place in relatively faint galaxies, too faint to be found in most flux limited surveys. 



In Figure El we show the preferred analytical LF fits of 
IBouwens et al.l (|2011cD and the data on which they are 
based, for samples at z ~ 5, 7 and 8. Their conclusion is 
that the characteristic knee in the luminosity function, 
L* , becomes slowly fainter with increasing redshift from 
4-8, and at the same time the faint end slope becomes 
steeper, reaching a value of a ~ 2 by z = 8. This is in ad- 
diti on to the overall nor malisation, </>*, lowering (but see 
also iMcLure et al1l2010l who find a fading L* and lower 
4>* but less evidence for a steepening slope). This indi- 
cates that an increasing proportion of star formation is 
occurring in fainter galaxies, and indeed formally the in- 
tegrated luminosity represented by a Schechter function 
with a > 2 diverges without some lower cut-off luminos- 
ity. However there are important caveats which pertain 
to these analyses: firstly, the LF parameters are based 
on entire samples, and any incompleteness or contami- 
nation (particularly difficult to rule out at the faint end) 
will introduce biases; secondly, even the HUDF is lim- 
ited to finding galaxies in the top few magnitudes of the 
LF, so that the conclusions about total star-formation 
rate (and hence the production rate of ionizing photons) 
are sensitive to the untested assumption that a Schechter 
function is the appropriate form, and to the large uncer- 
tainties on the measurement of the faint-end slope. 

We note that at high redshift the UV luminosity func- 



tion has limitations as a way of representing the whole 
population of galaxies, since they are typically only visi- 
ble when in the starbursting phase, and likely remain in 
this state for only a short duration based on the typi- 
cal ages, and apparent availability of molecular gas for 
star formation within them. However, particularly for 
understanding the contribution of galaxies to the reion- 
ization of the Universe, the measured UV luminosities, 
which are more representative of star formation occurring 
with the last ~ 15 — 100 Myrs, are the relevant quanti- 
ties. Furthermore, since the GRB itself is a sign of on- 
going massive star formation within the host, one would 
expect the likelihood that a GRB occurs to scale with 
the UV luminosity. Indeed, we would expect the region 
immediately underneath the burst to be extremely UV- 
bright, as is observed in the more local GRB population 
(|Fruchter et al.ll2006D . Hence we can meaningfully com- 
pare our limits on host emission with the observed galaxy 
LFs at z ~ 7. 

The 2cr luminosity limits for the GRB hosts are also 
shown graphically in relation to the high-z galaxy LFs 
in Figure [31 while their tabulated UV absolute magni- 
tude limits are given in Table [3J As can be seen, all the 
hosts are apparently below M* (the AB magnitude cor- 
responding to L*) at their respective redshifts, and in the 
case of our observations of GRB 090423 have the ability 



to probe to fainter limits than has so far been possible 
with the Ultra-Deep Field observations. This is because, 
despite the shallower depth of our images, we have prior 
knowledge that an object exists at this location, and so 
can accept a lower formal significance level for detection 
(since the chances of a random noise fluctuation are lower 
due to the much smaller area under consideration). Fur- 
thermore, since we do not need a blue "veto" filter we 
make more efficient use of the available exposure time; 
in contrast to some Lyman break searches which are lim- 
ited by the depth of their short- wavelength (blueward of 
Lya) imaging. 

These observations of GRB host fields allow us to test 
the validity of the galaxy LFs that have been derived 
from deep field observations. We assume that there is no 
dependence of GRB rate or luminosity on environmental 
parameters such as metallicity, which, as discussed in 
Section[TJ is certainly plausible at high redshifts when the 
bulk of star-forming galaxies were generally not highly 
enriched. We also assume that any dust extinction is 
minor, as argued in Section 13.1.11 Then the probability 
that a galaxy produces a GRB in some unit time (the 
probability of more than one GRB in realistic observing 
times being negligible) is approximately proportional to 
its rest-frame UV flux, since massive stars are responsible 
for producing both: i.e., Pgrb oc SFR oc Ltjv- Thus, 
because we would expect GRBs to be drawn randomly 
from the total stellar UV luminosity, this means that 
the host galaxies should be drawn from the luminosity- 
weighted galaxy LF. 

In practice, the best fit values for the faint end slope a 
and characteristic luminosity L* at each redshift consid- 
ered were determined f rom the fitting formulae given by 
iBouwens et al.l (|2011b[ ) , and are summarised in Table [31 
We can then calculate a luminosity-weighted luminos- 
ity function, which (suitably normalised) we take to be 
equivalent to a probability density function (PDF) for 
the intrinsic host luminosity: 



y(L) 



Lcj) z {L) 



(3) 



where 4> Z {L) is the luminosity function at redshift z. We 
emphasize again that, since we are going to compare to 
the fluxes measured in apertures at the exact GRB loca- 
tion, we are also assuming there is no significant offset 
between the GRB and its parent galaxy. 

Hence, the probability of observing a GRB in a galaxy 
of luminosity L < Lh os t can be obtained via the cumula- 
tive probability density function (CDF): 



Y(L 



J^ host L'<p z (L')dL' 



Setting a lower limit, L m - m , for the integral is phys- 
ically motivated since small dark halos (< 10 8 Af©) at 
z ~ 1 are expected t o retain little gas and form few 
stars (Re ad et al.l 12006). Given the apparent steepness 
of the faint end slope of the LF at z ~ 8 it is also im- 
portant that the value of -L m i n be sensibly chosen. It is 
not feasible to directly measure this lower limit at high 
redshift since it is well below the detection threshold for 
deep imaging, and will remain so even in the era of the 



James Webb Space Telescope (JWST) or ground-based 
Extremely Large Telescopes (ELTs). Howev er, it is pos- 
sible to provide estimates via simulations (Re ad et al.l 
2006), or from the star formation histories of the low- 
est m ass galaxies in the local universe (e.g. IWeisz et al.l 
2011) . This latter approach suggests that the lowest 
mass galaxies attain total stellar masses of ~ 10 6 M Q , 
over a Hubble time, implying mean star formation rates 
of ~ 10~ 5 — 10~ 4 M Q yr _1 . However, star formation is 
likely to be episodic, and characterized by periods when 
the SFR is markedly higher than this average, and other 
times when the star formation is inactive, and the galaxy 
effectively invisib le in the UV. Indeed, the study of 
IWeisz et all f|201 lh suggests that the SFR of local dwarfs 
was somewhat higher in the early Universe than their 
average SFR over the Hubble time. On the basis of this 
we adopt a cut-off value of L m - in = 4 x 10 23 ergs _1 Hz _1 
which is equivalent to Mab = —10 , and is similar to that 
considered by other recent studies (jBouwens et al.ll2011bl ; 
IKuhlen fc Fauche~ Gigucrc 2012). This corresponds to a 
star formation rate of ~ 5 x 10~ 4 M Q yr _1 . We note, 
that when a < 2 (i.e., for redshifts less than z ~ 7) the 
precise choice of L m in has little impact on the results. 

3.2.1. Analysis I 

We now consider the question: what is the probability 
that each host individually is fainter than the 2a upper 
bound we have inferred for its luminosity? To this end, 
we calculate Pia for each host by setting -Lhost m Equa- 
tion [4] equal to this upper bound. 

The results are summarised in Tableland illustrated 
in Figure dj which shows the CDFs for luminosity (ex- 
pressed relative to L* at the redshift in question) of two 
bursts from our sample, together with the 2a detection 
limits for the corresponding HST frames. We see that, 
given the above assumptions, it is not surprising to find 
individual hosts to be undetected since the majority of 
the likelihood lies below these 2a bounds. However, the 
joint probability that none of the hosts is detected is only 
P(none) = Y[ P2a — 0.17. While not a highly significant 
result, this does suggest that a non-evolving LF, in which 
more star formation was taking place in galaxies with lu- 
minosities around L* , would not sit comfortably with the 
apparent faintness of the GRB hosts. 

3.2.2. Analysis II 

A weakness of the above approach is that it does not 
make full use of the joint probability of the formally mea- 
sured fluxes at the positions of the whole sample of GRB 
hosts. As an alternative, we perform the following anal- 
ysis, again with the aim of test i ng the evolving galaxy 
LFs proposed by IBouwens et all (|2011bf ). 

As before, we construct a CDF for host luminosity, 
but now turn this into an equivalent CDF for observed 
flux density F using the cosmological luminosity distance 
for the given redshift, having accounted for the aperture 
corrections. The next step is to convolve this with the 
observational errors appropriate for the given GRB field 
observation, in order to obtain a CDF for the observed 
host galaxy flux density: 



Y(F) 



/~ G * (F'cj) z (F'))dF' 



(5) 
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Fig. 3. — Our 2<r upper limits on the magnitudes our indi vidual host galaxies (ve rtical dashed lines), compared to the luminosity functions 
of galaxies at z = 5, 7, 8 (solid curves), as determined by Bouwcns ct al. (2011c). In particular for GRB 090423, the power of deep HST 
observations to probe sources at known locations and redshifts to well below M* (indicated by upward arrows) can be readily appreciated. 
Note that GRB 090429B is plotted both in the top, highest redshift, panel appropriate for its best-fit photometric redshift, and also (lighter 
shade) in the middle, intermediate redshift, panel appropriate for the lower-limit of photo-z ss 6.5 (see text for more details). 



where G is a gaussian with a width a dictated by the sky 
noise in each image, measured from numerous sky aper- 
tures of equal radius to the source aperture. -Fobs is the 
formal flux measurement at the location of the GRB, and 
the minimum flux density, -F m i n , is appropriately scaled 
from the minimum luminosity, discussed above. 

For illustration, the CDFs for two of our bursts are 
shown in Figure [5] (red curves), along with the CDFs for 
true host flux for comparison (blue curves). The green 
lines indicate the formal measured flux density at each 
GRB position. Note, that for GRB 090423 we take a 



weighted average of the results for both filters, since they 
straddle 1500 A, while in the analysis of GRB 060927 we 
use the F110W flux density, being the closest match to 
1500 A. If the assumptions we have made are correct, 
then we would expect these measured flux densities to 
be drawn randomly and uniformly from to 1 on the 
cumulative probability axes. We can quantify this by 
calculating the average value for (Y(F b s )) = 0.46, con- 
sistent with an expected value of 0.5 ± 0.12 from the 
Central Limit Theorem. 
However, we can also repeat the analysis, but this time 
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Fig. 4. — These panels show the cumulative probability distribution s for the luminosities of two of our GRB h osts, based on assumin g 
GRB likelihood is proportional to UV luminosity, as described in Section l3.2l Galaxy LF parameters are taken from[Bouwcns ct al. ( 201 lbT l, 
The vertical lines correspond to the 2a detection limits for our HST images. As noted, in no case do we detect a host at this level of 
confidence. 
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Fig. 5. — Two examples of the cumulative probability density functi ons for the GRB host true flux densities (dashed curves; blue on-line), 
based on the luminosity- weighted galaxy LFs, as described in Section l3.2l These are then convolved with the observed measurement errors 
for each field to produce the solid (red on-line) curves which are the predicted CDFs for the measured flux densities. The dotted (green 
on-line) lines show the formal measured flux density at the location of the GRB in each case. Reading these across on the CDF axes we 
would expect to find locations drawn randomly from a uniform distribution between and 1. 



fix the LF parameters for all the hosts to be that found at 
lower redshift. Specifically, we first cho ose a = 1.73 and 
M* = -20.97, as measured at z w 3 bv lReddv fc Steidell 
(200<J. Thus we are testing here whether the hosts could 
be drawn from a luminosity function whose shape (but 
not normalisation) does not evolve from z ~ 3 to high 
redshift. The results in this case for the same two bursts 
are shown in Figure [6j Now all the measured flux den- 
sities at the GRB locations are close to or below the 
Y = 0.5 level, with a mean (Y(F \ 1S )) = 0.29, thus re- 
jecting the model at ss 96% confidence (~ 98% if we 



took z = 6.5 for the redshift of GRB 090429B and aver- 
aged the F105W and F160W limits). An alternative test 
wo uld be to fix th e para meters to those found at z = 6-7 
bv IMcLure eTaTl ((20l"l . a = 1.71 and M* = -20.08. 
Again, this model is weakly rejected at the « 90% level 
(«94%). 

Our result is not a highly significant, but does sup- 
port an evolving galaxy LF, with an increasing propor- 
tion of star formation occurring in faint galaxies. It also 
demonstrates that a larger sample and/or deeper limits 
on host emission can begin to provide important tests of 
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the high-z galaxy LF. 

3.3. Implications for reionization 

Whether UV radiation from star-forming galaxies is 
sufficient to bring about and sustain the reionization of 
the IGM a bove z ~ 6 is a long-standing question (e.g. 
lLoebl [20091 ). Recently a number of authors have argued 
that an increasingly steep faint end to the galaxy LF, 
motivated by theoretical considerations, may be able to 
achieve this, without resorting to extremely high Lyman 
continuum escape fractions (e.g. iBouwens et aLl|201lbt 
lLorenzoni et al.ll20lH Ikuhlen fc Faucher-Giguerell2012h . 
However, the form, steepness and faint-end cut-off of the 
galaxy LF at z > 7 are very poorly constrained by cur- 
rent data, since HST (and even in the longer-term JWST 
and ground-based 20-40 m class optical/nIR telescopes) 
can only directly probe the bright end of the LF. 

If GRBs are sampling star formation in an unbiased 
way, as we argue they may be at early times, then they 
provide an alternative window on the total star formation 
rate which would ultimately circumvent the necessity to 
detect the individual galaxies in which the star formation 
is occurring. Our results, even from the small sample 
of high-redshift bursts currently available, already sup- 
port an evolving galaxy LF over a non-evolving one, and 
therefore suggest that reionization may be brought about 
primarily from stars born in very faint proto-galaxies. 

3.4. Possible biases 

It is worth considering further possible physical effects 
which may be biassing our conclusions. If a significant 
amount of star formation at high-z is actually dust en- 
shrouded, which we believe is unlikely, then it would im- 
pact on the observability of GRBs and their hosts, as well 
as the LBG samples. One would generally expect that 
the GRBs for which afterglows are detected and redshifts 
estimated would typically be in the lower dust systems, 
if there is indeed a wide range, and so in that sense, their 
hosts could still be compared directly to the (also prefer- 
entially dust free) LBG samples. From the point of view 
of reionization, of course it is the low-dust star formation 
that is more likely to have a high escape fraction of ion- 
izing radiation, so our conclusions are likely to be valid 
in that respect. 

As discussed in section 1, we should also be concerned 
about possible GRB metallicity sensitivity. Both the- 
oretical considerations and observational evidence (the 
blue colours of the LBG samples) argue that few star- 
forming galaxies at z > 6 will have high metallicities 
(e.g. supersolar) and therefore the low rate of GRBs 
in such systems seen at low redshift is not likely to be 
an important factor at high-z. Of course, it could be 
that a contrary effect becomes important at some point, 
for instance if very low metallicity populations produce 
fewer and/or fainter GRBs. However, if that were the 
case then we would expect to lose, if anything, the lower 
mass halos, where any metals produced are most easily 
lost. Hence, such an effect would seem unlikely to result 
in finding an unusually faint population of hosts, and so 
again our basic conclusion would only be strengthened. 

4. CONCLUSIONS 

We have presented HST observations of a sample of 
six GRB host galaxies beyond z ~ 5. One host, that of 



GRB 060522 at z = 5.11 may be marginally detected, but 
the others are undetected to deep limits, typically Hab ~ 
28 at 2a. If GRBs are good tracers of the locations of star 
formation at high-z our results confirm that much, and 
probably the majority of, star formation then was taking 
place in small galaxies that are too faint to be detected 
even in the various HST deep field surveys. While the 
sample is small, the joint probability that none of the 
hosts are detected is consistent with a galaxy LF which 
is rapidly evolving to higher redshifts, and marginally 
inconsistent with an LF whose shape does not evolve. 

Our analysis does rely on two assumptions: that envi- 
ronmental conditions, such as abundance variations, do 
not produce appreciable variations in the SFR to GRB- 
rate ratio from galaxy to galaxy, and that dust content is 
generally negligible. Both of these are consistent with our 
current understanding of early galaxies, and of GRBs, 
but doubtless work is required to further clarify these is- 
sues. In terms of the effect of dust, we note that from 
the point of view of reionization, the un-corrected lumi- 
nosity function is really what we are interested in, since 
it is largely the unobscured star formation which will 
contribute to the inter-galactic UV radiation field. How- 
ever, if there were significant dust in some GRB hosts, 
it would weaken the connection between GRB-likelihood 
and apparent UV luminosity which we have assumed. 

This work demonstrates the potential power of GRB- 
selected galaxy samples to quantify the amount of star 
formation occurring in faint galaxies at early times, 
which is essential for understanding the budget of UV 
photons and their role in reionizing the Universe. In the 
future, deep imaging of larger samples of GRBs at high-z, 
combined with better understanding of any environmen- 
tal dependencies of GRB production, could provide much 
more stringent constraints on the faint-end of the galaxy 
luminosity function. Specifically, such a sample could in 
principle be used to fit for all the LF shape parameters, 
including the minimum cut-off luminosity, or, more em- 
pirically, simply determine the relative SFR fraction in 
bright (detected) galaxies compared to that in faint (un- 
detected) galaxies, without any prior assumption about 
the form of the LF. 

Just as this paper was submitted, t wo other papers on 
the same topic appeared as preprints. iBasa et all ()2012f ) 
have reported deep VLT observations of three high-z 
fields, of GRBs 060522, 060927 and 080913. The point- 
source magnitudes reached are rather shallower than the 
HST limits, and t heir non-dete c tions therefore consistent 
with our results. iTrenti et~ al. (2012) instead performed 
a theoretical analysis, aimed at predicting the fraction 
of GRB hosts expected to be detected in deep HST 
imaging at different redshifts, under various model as- 
sumptions. For example, they predict that 50% of GRB 
hosts at z ~ 5 should be detected in a survey reach- 
ing Mab ~ — 18, consistent with our possible detection 
of the GRB 060522 host. To the same rest-frame limit, 
they predict fewer than 20% of hosts should be detected 
at z ~ 8, again, in agreement with our findings. 

Traditional blank-field imaging surveys, while poten- 
tially finding larger samples of candidate galaxies, are 
limited to detecting just the relatively brightest exam- 
ples at the most extreme redshifts. This may be true 
even in the era of the JWST, especially given the ap- 
parent dearth of z ~ 10 candidates located so far in the 



12 



£T 6 



-1 1 1 


1 1 — 


— 1 1 1 1 


1 1 1 1 1 1— 


- 








- 




/ 








/ / 








/ / 








/ / 




- 




~ 1 / 


- 






■1 / 








\ 


050904 . 


— 1 1 1 


1 1 — 


— 1 1 1 1 






-50 



50 
F(n3y) 



100 -50 



50 
^(nJy) 



100 



Fig. 6. — As with Figure[5]but here we have fixed the shape of the input galaxy LFs to a = 1.73, M* = —20.97 which is representative 
of the population at 2 S3 3. The shallower slope and brighter characteristic luminosity of this LF leads to our faint flux limits corresponding 
to systematically lower values of the CDF, Y(F). 
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HUDF ()Oesch et al.ll2012D , which could be indicating an 
even more rapid evolution of the galaxy LF parameters. 
The approach we, and the other studies mentioned above, 
have adopted provides a crucial complementary insight 
into early galaxy evolution. 
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TABLE 2 

Log of HST observations of the host galaxies of GRBs at z > 5 



Date 


UT Time 


Filter 


\ 1 A "i 

Arest 


Exp (s) 


^obs (nJy) 


AB mag limit 


M A/(l+z) 


060522 
















17 Oct 2010 


10:30 


F110W 


1888 


8395 


7 ± 4 


> 28.13 


> -18.35 


060927 
















29 June 2007 


11:30 


F160W 


2376 


10240 


7 ± 5 


> 27.75 


> —18.84 


25 Sept 2010 


14:30 


F110W 


1782 


13992 


4±3 


> 28.57 


> -18.02 


050904 
















9fi CInnt 900^ 


91 -nt 




1 97Qt 


4216 


q _|_ got 




, iq nc: 


uouyio 
















30 Nov 2009 


16:10 


F160W 


1988 


7818 


3 ± 6 


> 27.92 


> -19.00 


0QD42S 
















24 Jan 2010 


11:34 


F160W 


1665 


13029 








25 Jan 2010* 


14:44 


F160W 


1665 


13029 


4±3 


> 28.36 


> -18.88 


26 Jan 2010* 


13:06 


F125W 


1353 


13029 








27 Jan 2010* 


13:04 


F125W 


1353 


13029 








22 Oct 2010 


18:23 


F125W 


1353 


13029 








27 Oct 2010 


16:36 


F125W 


1353 


13029 


-2 ±2 


> 30.29 


> -16.95 


090429B 
















10 Jan 2010 


21:54 


F160W 


1478 


2412 








22 Feb 2010 


19:22 


F160W 


1478 


2412 


7±5 


> 27.78 


> -19.65 








(2049) 








(> -19.09)* 


24 Feb 2010 


03:19 


F105W 


1014 


2412 








28 Feb 2010 


13:56 


F105W 


1014 


2412 


-1±5 


> 28.49 


> -18.73 








(1407) 








(> -18.17)* 



Note. — Log of HST observations of the locations of our sample of bursts giving the details of the 
observations and the pivot wavelength of the filter at the assumed redshift, A res t. Fobs represents the 
measured (sky-subtracted) residual flux density in an aperture of diameter f//4, centred on the position 
of the afterglow. The 2ct point-source limits on apparent magnitude (after aperture correction) and 
corresponding absolute magnitude for each filter are also given. *Images affected by persistence, t Filter 
central wavelength and flux density calibration have been corrected for the fact that the Lya break for 
GRB 050904 falls roughly in the middle of the filter response, and so in effect only the longer wavelength 
part of the filter passband is actually sensitive to the flux from the host, and the quoted limit is for the 
flux density redward of Lyc?. * Result if GRB 090429B was actually at the lower limit for its photometric 
redshift of z as 6.5. 



TABLE 3 

Assumed luminosity function parameters at each redshift, and derived host properties and 

probabilities 



Burst z M* VVAB a m* F*(njy) L host /L* SFR(Moyr -1 ) P 2ct Y(F oba ) 



060522 


5.11 


-20.59 


1.80 


26.22 


118 


< 0.09 


< 0.88 


0.55 


0.47 


060927 


5.47 


-20.49 


1.81 


26.43 


97 


< 0.08 


< 0.65 


0.55 


0.43 


050904 


6.29 


-20.26 


1.85 


26.83 


67 


< 0.59 


< 4.1 


0.94 


0.28 


080913 


6.73 


-20.14 


1.88 


27.33 


43 


< 0.21 


< 1.3 


0.82 


0.44 


090423 


8.23 


-19.72 


1.95 


27.97 


24 


< 0.09 


< 0.38 


0.78 


0.30 


090429B 


9.4 


-19.39 


2.01 


28.59 


13 


< 0.78 


< 2.4 


0.96 


0.82 




(6.5)* 


(-20.20) 


(1.87) 


(26.68) 


(77) 


(< 0.13) 


(< 0.84) 


(0.71) 


(0.46) 



Note. — Columns 3—6 give the galaxy luminosity function pa rameters assumed at the redshifts of each 
GRB obtained from the fitting formulae of Bouwcns ct al. (2011b), both in absolute terms and translated to 
apparent values (including correction to finite aperture) . Columns 7—8 list the 2a limits on host luminosity and 
star-formation rate; in the case of GRB 060927 based on just the F110W image, and for GRB 090423 based 
on the weighted average of both filters, as described in the text. Finally, columns 9-10 give the calculated 
probability that the host would be below this 2a limit, P^a and the position in the cumulative prob ability 
density function of the observed flux density, Y(-Fob s )j based on an evolving galaxy LF (Section 13.21 1. The 
combined probability of the 2a upper limits, the product of the figures in column 9, is n^o- =0.17. * Result 
if GRB090429B was actually at the lower limit for its photometric redshift of z 6.5. 
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